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Abstract 

We consider the polynuclear growth (PNG) model in 1 + 1 dimen- 
sion with flat initial condition and no extra constraints. Through the 
Robinson-Schensted-Knuth (RSK) construction, one obtains the mul- 
tilayer PNG model, which consists of a stack of non-intersecting lines, 
the top one being the PNG height. The statistics of the lines is trans- 
lation invariant and at a fixed position the lines define a point process. 
We prove that for large times the edge of this point process, suitably 
scaled, has a limit. This limit is a Pfaffian point process and identical 
to the one obtained from the edge scaling of Gaussian orthogonal en- 
semble (GOE) of random matrices. Our results give further insight to 
the universality structure within the KPZ class of 1 + 1 dimensional 
growth models. 

1 Introduction 

The polynuclear growth (PNG) model is the best-studied growth model from 
the KPZ class in one spatial dimension. Since in one dimension the dynamical 
scaling exponent is z = 3/2, the correlation length increases as t^/^ for large 
growth time t. The exponent is universal, but different classes of initial 
conditions lead to distinct scale invariant statistical properties of the surface 
in the large t limit. To go beyond the exponents and to determine the exact 
scaling functions one has to analyze some solvable models. In this paper 
we consider the PNG model. The surface height at time t is denoted by 
X I— > h{x, t). On the surface new islands of height one are created at random 
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with intensity q. The islands spread with unit speed and simply merge upon 
contact. 

For the PNG model the exact scaling function for stationary growth is 
known, see [Hj, relying on previous results by Baik and Rains jHEj- In the 
case where initially h{x, 0) = and nucleations are constrained to occur only 
above a first spreading layer, the surface has typically the shape of a droplet. 
In this geometry, for fixed but large t, the spatial statistics of the surface is 
well understood. Subtracting the deterministic part, it is proved in [12] that 
the self-similar shape fluctuations are governed by the Airy process. 

In numerical simulations, one starts the growth process mostly with a 
flat substrate, i.e., 0) = with no further constraints. Thus it would 
be of interest to understand the statistics of a; i— > h{x,t) at large t. It is a 
space translation invariant process and the only available result |31 is the 
one-point distribution, 

limP(/i(0,t) <2t + ^f/^) =^1(^2^/3), (1.1) 

with Fi the GOE Tracy- Widom distribution (here ^ = 2). The limiting 
function Fi is linked to the Gaussian Orthogonal Ensemble (GOE) of random 
matrix theory as follows. Let us denote by Xraa^,N the largest eigenvalue of a 
N X N GOE matrix, in units where E(Amax,Af) = 2iV. Then 

lim P (A^ax,iv <2N + = Fi(0, (1.2) 

see |2ni- -^1 is given in terms of a Fredholm determinant, compare ()2.1()|1 
below. A plot of dFi(^)/d^ in semi-logarithmic scale is available in [T^ . 

The result fll.l|) leaves open the joint distribution at two space points, 
even more the full process with respect to x. From the general KPZ scaling 
theory, see JHI for an exposition, a meaningful limit is expected only if the 
two points are separated by a distance of order t^/^. Thus the issue is to 
determine the limit 

lim P (h{0, t)<2t + ^if/^, h{Tt^/^, t) < 2t + ^2^^^^) = ? (1.3) 

Of course, the marginals are Fi (^i22/3) and Fi(^22^/^). But this leaves many 
choices for the joint distribution. 

In the present contribution we will not succeed in removing the question 
mark in p.3|) . However, we will make a big step towards a well-founded 
conjecture. The idea to progress in the direction constituting the main body 
of our paper was set forward by Kurt Johansson in a discussion taking place at 
the 2003 workshop on growth processes at the Newton Institute, Cambridge 
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(as communicated by Herbert Spohn). In a somewhat rough description, 
underlying the PNG process there is a hne ensemble constructed through 
the Robinson-Schensted-Knuth (RSK) algorithm. Its top line at time t is 
the height x ^— h{x,t). Thus the question mark in refers to the top 
line. But instead of (jl.Hj) we will study, for large t, the line ensemble at 
fixed X = close to the top line. As in a way already suggested by the Baik 
and Rains result, under suitable scaling the lines at x = have indeed the 
statistics of the top eigenvalues of GOE random matrices. The implications 
for will be discussed after explaining more precisely our main result. 

2 Main result 

First we describe the PNG model with flat initial conditions, and secondly 
recall some random matrix results on GOE eigenvalues, as needed to state 
the scaling limit. 

2.1 Polynuclear growth (PNG) model and RSK con- 
struction 

The polynuclear growth (PNG) model considered here is a 1 + 1 dimensional 
model. One way to view PNG is via a graphical construction involving 
Poisson points. Consider a fixed T > and let u; be a countable configuration 
of points in R X [0,T]. For any compact subset B of H x [0,T], denote the 
number of points of u; in i? by n{B){uj). Then 

Q = {uj\n{B){uj) < oo,V compact B cRx [0,T]} (2.1) 

is the set of all locally finite point configurations in R x [0, T]. The Poisson 
process with intensity ^) > in R x [0, T] is given by setting the probability 
Pt such that 

PT{{n{B) = n}) = ^^l^l)" e-gl^' (2.2) 

nl 

for all compact B C IRx [0, T], and the family of random variables {n{Bj),j = 
1, . . . , m} with Bi nBj = ^ for i 7^ j, is always independent. In what follows 
we set g = 2. 

For each u E flwe define the height function h{x, t){u!), {x, t) G Rx [0, T], 
by the following graphical construction. Because of fiat initial conditions, we 
set h{x,0){u!) = and we call nucleation events the points of u. Each 
nucleation event generates two lines, with slope +1 and —1 along its forward 
light cone. A line ends upon crossing another line. In Figure ^ the dots 
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t = T 



Figure 1 : Graphical construction generating the surface height from the Pois- 
son points. 



are the nucleation events and the hues follow the forward light cones. The 
height h{x,t){uj) is then the number of lines crossed along the straight path 
from (x, 0) to (x, t). Since u is locally finite, it follows that x h{x,t){uj), 
t G [0,T], is locally bounded and the number of discontinuities is locally 
finite. 

The interpretation of the graphical construction in terms of a growing 
surface is the following. The surface height at position x G R and time t > 
is h{x, t) e The initial condition is h{x, 0) = for all x G R. For fixed time 
t, consider the height function x ^ h{x,t). We say that there is an up-step 
(of height one) at x if h{x, t) = lim^^i^. h{y, t) + 1 and a down-step (of height 
one) at x if h{x,t) = liiriyix h{y,t) + 1. A nucleation event which occurs at 
position X and time t is a creation of a pair of up- and down-step at x at time 
t. The up-steps move to the left with unit speed and the down-steps to the 
right with unit speed. When a pair of up- and down-step meet, they simply 
merge. In Figure^the dots are the nucleation events, the lines with slope — 1 
(resp. -|-1) are the positions of the up-steps (resp. down-steps). Other initial 
conditions and geometries can be treated in a similar fashion. For example, 
if h{x, 0) is not for all x, it is enough to add additional lines starting from 
the t = axis with slope ±1 reflecting the up/down direction of the steps 
at t = 0. Another interesting situation is the PNG droplet, where one starts 
with fiat initial conditions and there are no nucleation points outside the 
forward light cone starting at (0,0). 

To study the surface height at time T, x h{x,T), it is convenient to 
extend to a multilayer model. This is achieved using the RSK construction. 



4 




We construct a set of height functions h^{x, t){uj), {x,t) G IR x [0, T], £ < as 
follows. At t = we set hi{x, 0) = i with i = 0, —1, . . ., £ denoting the level's 
height. The first height is defined by ho{x,t){uj) = h{x,t){ijj). The meeting 
points of the forward light cones generated by the points of uo are called 
the annihilation events of level 0. h_i[x,t){oj) is constructed as hQ{x,t){(jj) 
but the nucleation events for level —1 are the annihilation events of level 
and h_i{x,t){uj) + 1 equals the number of lines for level —1 crossed from 
(x, 0) to (x, t). In Figure 121 the nucleation events of level —1 are the empty 
dots, whose forward light cones are the dotted lines. Setting the annihilation 
events of level j as the nucleation events for level j — 1, the set of height 
functions he{x,t){uj) is defined for all i < 0. The line ensemble for t = T, 
i.e., {hi{x,T), i < 0} is represented in FigureEl 

The point process which describes this line ensemble at z = is denoted 
by (^'^^ and given by 

Aflat/ N ^ f 1 if a line passes at (0, j), ^2 3) 

1^ if no line passes at (0,j). 

From the Baik and Rains result we know that the largest j such that Ct^*^) 7^ 
is located near 2T and fluctuates on a T^/^ scale. The edge rescaled point 
process is defined as follows. For any smooth test function / of compact 
support 

V^if) = J2^T\j)f{ij - 2T)/(TV32-2/3)), (2.4) 
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Figure 3: Line ensemble for t = T for the point configuration of Figure |21 



the factor 2~^/'^ is the same as in (jl.lj) . Notice that in ()2.4|1 there is no 
prefactor to the sum. The reason is that close to 2T, the points of (^'^^ are 
order T^/^ apart and 77^*^* remains a point process in the limit T ^ 00. r^f.^* 
has a last particle, i.e., rj^^^i^) = for all ^ large enough, and even in the 
T — ^ cxD limit has a finite density which increases as y/—^ as ^ —>■ —00. 
Consequently the sum in ()2.4|1 is effectively finite. 

2.2 Random matrices 

The Gaussian Orthogonal Ensemble (GOE) is the set of N x N real 
symmetric matrices distributed according to the probability measure 
Z^^exp[-TrMV2A^]dM, with dM = Ui<i<j<N The eigenvalues are 

then distributed with density 

1 ^ 
^ n |A.-A,inexp[-A^/2iV]. (2.5) 

^ ^<i<j<N i=l 

Let us denote by C_^*^^ the point process of GOE eigenvalues, i.e., C^*^^(a;) = 
SjLi ^{^ ~ Aj). At the edge of the spectrum, 2N, the eigenvalues are order 
N^^^ apart. The edge rescaled point process is then given by 

r^GOE(^) ^ iVV3^GOE(2jv + (2.6) 
and for / a test function of compact support. 



N 

if) = - 2iV)/iVV3) = / dUXOvT^iO- (2.7) 



We denote by rj'^'^^ the limit of rj'^'^^ as N 00 
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The limit point process rj'^'^^ is characterized by its correlation functions 
as follows. Let us denote by Pgoe(^1' • • • 5 Cn) the n-point correlation functions 
of T]'^'^^, i.e., the joint density of having eigenvalues at ^1, . . . , ^n- Then 

pSeI^I, . . . , en) = Pf 0)].,=l,...,n (2.8) 

where Pf is the Pfaffian and G"^*-*^ is a 2 x 2 matrix kernel with elements 

POO 

G^:?^(6,6) = / dAAi(a + A)Ai'(6 + A)-(a--6), (2.9) 

poo 1 POO 

^^^"^(^1,6) = J dAAi(ei + A)Ai(6 + A) + -Ai(6)y dAAi(6-A), 
G??^(ei,6) = -G^^^(6,ei) 

1 POO poo 

Gg°''(6,6) = 4 dA d/iAi(6-/i)Ai(ei-A)-(6-6), 

and Ai is the Airy function IJ. The notation (^i ^ ^2) nieans that the 
previous term is repeated with and ^2 interchanged. For an antisymmetric 
matrix A, Pf (A) = ^Det(v4), see (glU) for the definition of the Pfaffian. The 
GOE kernel was studied in jSHI- It is not uniquely defined, for example the 
one reported in jHJ CB] differs slightly from the one written here, but they 
are equivalent because they yield the same point process. The point process 
^GOE jg uniquely determined by its correlation functions [2U]. 

Finally let us remark that Fi can be written in terms of a Fredholm 
determinant 

F,{0 = Jhn E^(n(l-l(C,oo)((A,-2Ar)/iVV3))'j 

= VDet(l- J-iGGOE)^ (2.10) 

where J is the matrix kernel J = ^ ) ' "^^^ determinant in ()2.1()j) is 

the Fredholm determinant of the kernel K = J^^G'^'^^ on the measure space 
{{^,00) X {1,2}, da; x z/) with dx the Lebesgue measure and z/ the counting 
measure on {1,2}, i.e., 

^etil-K) = J2^—^ dei---de„ Det[ir.„,(e.,6)].,=i,...,„- 
k=0 -^K'-)" n,...,i„e{i,2} 

(2.11) 

Remark: One can also consider instead of Det(l — K) the determinant 
Det(l — K) with K the operator with kernel K. K is not trace-class on 



7 



L^((^, oo), dx) © L'^iiC,, oo), dx) because it is not even Hilbert-Schmidt. Nev- 
ertheless it is possible to make sense of it as follows. K is Hilbert-Schmidt 
in the space oo), 6'~^dx) © oo), 6'dx) where 9 is any positive 

weight function growing at most polynomially at x — > oo and satisfying 
e L^{{^, oo), dx). Moreover the sum of the diagonal terms, Tt{K), is ab- 
solutely integrable. Then the modified Fredholm determinant, which has also 
the series development (jTTT|l . is defined by Det(l-it) = e''^^^) Det2(l-it) 
with Det2 the regularized determinant [Zj. This functional analysis point of 
view is used by Tracy and Widom [22] to show the convergence of the mod- 
ified Fredholm determinants in the N —>■ oo limit. 

2.3 Scaling limit 

As our main result we prove that the point process r]^^^ converges weakly to 
the point process rj'^'^^ as T oo. 

Theorem 2.1. For any m G N and smooth test functions of compact support 

fit ■ ■ ■ 1 fni! 

^lim ET(nr^r(^)) =E(nr^««^(/.)). (2.12) 

refers to expectation with respect to the Poisson process measure Pt- 
The expected value on the r.h.s. of ()2.12|1 is computed via the correlation 
functions (Q. 

As announced in the Introduction, the result of Theorem 12.11 is a first 
step towards a conjecture on the self-similar statistics of the PNG with flat 
initial conditions. The starting observation is that, as for the PNG, also to 
random matrices one can introduce a line ensemble in a natural way. Let M 
he a N X N random matrix in the GOE, resp. GUE, ensemble. As noticed 
by Dyson [S] when the coefficients of M are independent Ornstein-Uhlenbeck 
processes, then the eigenvalues Xj{t) of M = M{t) satisfy the set of stochastic 
differential equations 

dA,(<) = ( - ^X,it) + f ± JJfi^) * + . . = 1. .... N, 

(2.13) 

with {bj{t), j = 1, A^} a collection of A^ independent standard Brownian 
motions, /? = 1 for GOE and (3 = 2 for GUE. We refer to the stationary 
process of ()2.13j) as Dyson's Brownian motion. Note that for 13 >1 there is 
no crossing of the eigenvalues, as proved by Rogers and Shi |17j . 
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Let us denote by the point process of GUE eigenvalues, i.e., 

(^^^{x,t) = J2f=i ^(^ ~ ^ji't)) where Xj{t) are the GUE eigenvalues at time 
t. At the edge of the spectrum the eigenvalues are order N^^^ apart and the 
space-time edge rescaled point process is given by 

r^GUE^^^ ^) ^ iVi/3cGUE^2iV + ^N'/^, 2riV2/3). (2.14) 

Its limit as — * CX3 is denoted by rj'^^^iC,, t). r^*^UE jg ^ determinantal point 
process and its space-time kernel is the extended Airy kernel. The top line of 
^GUE -g gjygj-^ ^]^g Airy process, denoted by A{t), which appears in [T5| l9] 
with more detailed properties investigated in [21 ESI 1^ l2Zj- The height 
statistics x h-* h{x, t) for the PNG droplet is linked to the Airy process by 

hm T-^/^(h{TT^/\T) - 2^T2 - {rT'^/^y^ = A{t), (2.15) 

where the term subtracted from h is the asymptotic shape of the droplet J3] . 
To obtain this result, Prahofer and Spohn consider the line ensemble obtained 
by RSK and define a point process like ()2.3p but extended to space-time. It is 
a determinantal point process and in the edge scaling it converges, as T ^ oo, 
to the point process associated with the extended Airy kernel. Thus they 
prove not only that the top line converges to the Airy process, but also that 
the top lines converge to the top lines of Dyson's Brownian motion with 
(3 = 2. 

One can extend r/^*^^ of ()2.6|1 to space-time as in ()2.14|1 . The conjecture 
is that, under edge scaling, the process x i— >• h{x,T) for fiat PNG is in 
distribution identical to the largest eigenvalue of Dyson's Brownian motion 
with (3 = 1. The result of Theorem 12.11 makes this conjecture more plausible. 
In fact we now know that, not only /i(0, T) in the limit T ^ oo and properly 
rescaled is GOE Tracy- Widom distributed, but also that the complete point 
process rj^^ converges to the edge scaling of Dyson's Brownian motion with 
(3 = 1 for fixed time. For (3 = 1 Dyson's Brownian motion one expect that 
under edge scaling the full stochastic process has a limit. More explicitly, one 
focuses at the space-time point (2A^, 0), rescales space by a factor N'^^'^, time 
by A^2/3^ g^j^jj expects that the statistics of the lines has a limit for ^ oo. 
It could be that this limit is again Pfaffian with suitably extended kernel. 
But even for (3 = 1 Dyson's Brownian motion this structure has not been 
unravelled. 

The rest of the paper is organized as follows. In Section 01 we explain 
the line ensemble which will be used to obtain our result. It differs from the 
one of Figure El The end points of the line ensemble gives a point process, 
whose correlation functions are obtained in Section |31 They are given in 
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Figure 4: A configuration with three Poisson points in the triangle A+ and 
their sjTiimetric images with respect to the t = axis. The path 7 is the 
bold line. On the right we draw the top lines of the line ensemble associated 
to 7, {H,{s),j<0,se[0,TV2]}. 



term of a 2 x 2 matrix kernel which is computed for fixed T in Section El 
Section IHl is devoted to the edge scaling of the kernel and its asymptotics. 
Finally Section [7| contains the proof of Theorem 12.11 



3 Line ensemble 

3.1 Line ensemble for the S symmetry 

The line ensemble for flat PNG generated by RSK at time t = T is not easy 
to analyze because there are non-local constraints on the line configurations. 
Instead, we start considering the point process C^^^. First remark that this 
point process depends only on the points in the triangle A+ = {{x,t) G 
R X ]R+|t e [0,T], \x\ < T — t}. We then consider the Poisson points only 
in A+ and add their symmetric images with respect to the t = axis, which 
are in A_ = {{x,t) E R x R^\{x, -t) e A+}. We denote by Ct"" the point 
process at (0, T) obtained by RSK construction using the Poisson points 
and their symmetric images, see Figure El To study we consider a 
different line ensemble. Let us consider the path in space-time defined by 
7(5) = (T — s/\/2, s/\/2), s G [0,Ta/2], and construct the line ensemble 
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{Hj(s),j <0,sE [0,T^/2]}, as follows. The initial conditions are Hj(0) = j 
since the height at t = is zero everywhere. Every times that 7 crosses a 
RSK line corresponding to a nucleation event of level j, Hj has an up-jump. 
Then the point process Ct^ given by the points {Hj(T\/2), j < 0}. In 
Proposition 13.21 we show that (^^^ can be recovered by Ct'™, in fact we prove 
that hj{0,T) = 1{Hj{TV2)+j). 

Next we have to determine the allowed line configurations and their dis- 
tribution induced by the Poisson points. This is obtained as follows. We 
prove that the particle-hole transformation on the line ensemble {Hj{s),j < 
0, s e [0,Ta/2]} is equivalent to a particular change of symmetry in the po- 
sition of the nucleation events, and we connect with the half-droplet PNG 
problem studied by Sasamoto and Imamura [THl . 

Young tableaux 

Let cr = ((t(1), . . . , a{2N)) be a permutation of {1, . . . , 2N} which indicate 
the order in which the Poisson points are placed in the diamond A+ U A_. 
More precisely, let {xi,ti) be the positions of the points with the index i = 
1, . . . ,N such that -|- Xi is increasing with i, and a is the permutation such 
that to-(j) — is increasing in i too. Let us construct the line ensembles 
along the paths (T, 0) (0,T) and (— T, 0) (0,T). The relative position 
of the steps on the line ensembles are encoded in the Young tableaux S'(cr) 
and T(cr) constructed using Schensted's algorithm. If the fc*^ step occurs in 
line Hj, then in the Young tableau there is a in row j, see Figure El 

In our case the points are symmetric with respect to the axis t = 
and we refer to it as the symmetry S. In the case studied in ^], the 
points are symmetric with respect to the axis x = and we call it the 
symmetry 0. Consider a configuration of points with symmetry S and let a 
be the corresponding permutation. The RSK construction leads to the line 
ensembles of T{a) and S{a) as shown in the left part of Figure |31 If we apply 
the axis symmetry with respect to x+t = 0, then we obtain a configuration of 
points shown in the right part of Figured The points have now the symmetry 
IZl and the corresponding permutation a is obtained simply by reversing the 
order of a, that is, if a = . . . , (r{2N)) then a{j) = a{2N + 1 - j)- By 

Schensted's theorem jTH], 

S{a) = S{a)\. (3.1) 

Moreover, the positions of the steps in the line ensembles of S{a) and S{(y) 
occurs at the same positions, but of course in different line levels. Figure El 
shows an example with cr = (2 4 5 1 6 3), for which the Young tableaux are 
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Figure 5: Line ensembles for a = (24516 3) and a = (3 61542). For 
the even levels, i = 0, —2, we use the solid lines and for the odd levels, 
i = —1, —3, the dashed lines. The line ensemble of S{a) corresponds to the 
line ensemble {Hj{s)J <0,s e [0, Tv^]} of Figure H 



Particle-hole transformation 

At the level of line ensemble we can apply the particle-hole transformation, 
which means that a configuration of lines is replaced by the one with jumps at 
the same positions and the horizontal lines occupy the previous empty spaces, 
as shown in Figure IHl Let us start with the line ensemble corresponding to 
S{a), then the Young tableau for the hole line ensemble is given by 5'(o")*. In 
fact, the information encoded in S{a) tell us that the j^^ particle has jumps 
at (relative) position S{a)j^k, > 1, for j > 1. On the other hand, the 
jth i^qIq j^g^g jumps where the particles have their j^^ jump. Therefore the 
particle-hole transformation is equivalent to the symmetry transformation 
S ^ 0. 

Allowed line configurations and measure 

Sasamoto and Imamura pH] study the half-droplet geometry for PNG, where 
nucleation events occurs symmetrically with respect to a; = 0, i.e., with the 
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Figure 6: Particle (solid) and hole (dashed) line ensembles for the example 
of Figure El The particle line ensemble is the one associated with S{a), and 
the hole line ensemble is the one of S{a) reflected with respect to the line 
j = 1/2. The pairing rule is shown by the brackets. 



symmetry. In particular, they prove that the point process at x = 
converges to the point process of eigenvalues of the Gaussian Symplectic 
Ensemble (GSE). Its correlation functions have the same Pfaffian structure 
as GOE but with a different kernel. In a way the line ensemble they study is 
the hole line ensemble described above, thus their edge scaling focuses at the 
top holes, i.e., in the region where the lowest particles are excited. Notice 
that the change of focus between particles and holes changes the statistics 
from GSE to GOE. This differs from the case of the PNG droplet ^H] where 
for both holes and particles the edge statistics is GUE. Although the result 
of JHl cannot be applied directly to our symmetry, some properties derived 
there will be of use. 

From we know that for the symmetry a hole line configuration is 
allowed if: a) the lines do not intersect, b) have only down-jumps, c) they 
satisfy the pairing rule: H^f^{Ty/2) = H^°\{T^) for all j > 1. This 
implies that for the symmetry S a line configuration {Hj} is allowed if: a) 
the lines do not intersect, b) have only up-jumps, c) Hj{T\/2) — Hj{0) is 
even for each j < 0. Moreover, there is a one-to-one correspondence between 
allowed configurations and nucleation events. The probability measure for 
the line ensemble turns out to have a simple structure. Consider Poisson 
points with intensity g and symmetry S. Each Poisson point (x, t) G A+ has 
a probability gdxdt of being in [x, x+dx] x [t, t+dt]. In the corresponding line 
ensemble this weight is carried by two jumps, therefore the measure induced 
by the points on a line configuration {Hj} is given by ^/g^-^^^^^ ™ ^^^^ times 
the uniform measure. 
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3.2 Flat PNG and line ensemble for S symmetry 

The correspondence between the point process (^"^^ and Ct'™ is as fol- 
lows. Let us consider a permutation a with Young tableau S{a) of shape 
(Ai, A2, . . . , Am)- Let, for k < m, 0^(0") be the length of the longest subse- 
quence consisting of k disjoint increasing subsequences. 

Theorem 3.1 (Greene [8J). For all k = 1, . . . ,m, 

ak{a) = Xi + --- + Xk. (3.2) 

The geometric interpretation is the following. Let a be the permutation 
which corresponds to some configuration of Poisson points in A_|_ U A_. Then 
Ofc is the maximal sum of the lengths of k non-intersecting (without common 
points) directed polymers from (0, — T) to (0,T). 

Proposition 3.2. Let tx he a Poisson point configuration in A_(_ and let the 
corresponding Young tableau S{ti) have shape (Ai, A2, . . . , A^). Let n be the 
configuration of points on U A_ with symmetry S which is identical to it 
in A-|.. Then S{n) has shape (Ai, A2, . . . , Am) = (2Ai, 2A2, . . . , 2Am)- 

Proof. To prove the proposition is enough to prove that cifc(vr) = 2afc(7r) for 
k = 1, . . . ,m. 

i) ak{fc) > 2afc(7r): it is obvious since we can choose the k directed polymers 
on TT by completing the ones on vr by symmetry. 

ii) afc(7f) < 2ak{7i): assume it to be false. Then there exists k directed 
polymers in A_|_ and k in A_ such that the total length is strictly greater 
than 2afc(7r). This implies that at least one (by symmetry both) of the sets 
of k directed polymers has total length strictly greater that afc(7r). But this 
is in contradiction with the definition of akiir), therefore afc('7r) < 2afc(7r). □ 

Since Ai_j = hj(0,T) — j and Ai_j = Hj{T^/2) — j, it follows from this 
proposition that 

h,iO,T) = liH,iTV2) + j) (3.3) 

for all j < 0. 

4 Correlation functions 

Non-intersecting lines can be viewed as trajectories of fermions in discrete 
space Z and continuous time [0,T-\/2]. Let us start with a finite number of 
fermions, 2N, which implies that only the information in the first 2N levels 
in the RSK construction is retained. For any configuration, the number 
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of non perfectly fiat lines, is obviously bounded by the number of Poisson 
points in A_(_. On the other hand for fixed T, the probability of having a 
number of Poisson points greater than 2N decreases exponentially fast for 
large. First we derive an exact formula for the n-point correlation function 
for finite A^, and then take the limit iV — >■ cx) so that, for any fixed T, each 
line configuration contains all the information of the Poisson points. Finally 
we consider the asymptotic for large T. 

4.1 Correlation functions and Pfaffians 

The correlation functions of the point process Ct^™ turn out to be expressed 
as Pfaffians. Therefore we first review the Pfaffian ensemble. Let A = 
[Aj jjj . be an antisymmetric matrix, then its Pfaffian is defined by 

N 

pf(^)= E (-i)''^'n^— - (4.1) 
CTGS2JV «=i 

<^2i-l<0'2i 

where S2N is the set of all permutations of {1, ... , 2N}. Notice that the Pfaf- 
fian depends only on the upper triangular part of A. For an antisymmetric 
matrix the identity Pf (A)^ = Det(A) holds. 

The Pfaffian ensemble is introduced in jTHI, see also j5T]. Let {X,fi) be 
a measure space, /i, . . . , f2N complex- valued functions on X, e{x, y) be an 
antisymmetric kernel, and define by 

. . .,X2n) = -77—^(i^[fj{Xk)]j,k=l,...,2NPi[£{Xj,Xk)]j,k=l,...,2N (4.2) 
^2N 

the density of a 2iV-dimensional probability distribution on X"^^ with respect 
to fi'^'^^ , the product measure generated by fi. The normalization constant 
is given by 

Z2N = [ d2%Det[/,(xfc)],,fc=i,..„2^Pf[£(x„a;fc)],,,=i,...,2N = (2iV)!Pf[M] 

(4.3) 

where the matrix M = [Mjj]jj=i^ ^27v is defined by 

Mi,j= fi{x)e{x,y)fj{y)dfi{x)dfi{y). (4.4) 

The n-point correlation functions p'^"^(xi, . . . , Xn) of a point process with 
measure ()4.2j) are given by Pfaffians 

p^"\xu . . . ,x„) = Ff[K{xi,Xj)]ij=i,,..,n (4.5) 
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where K{x,y) is the antisymmetric kernel 



K{x,y)={ ^''f^y\ ^''f^y\ ] (4.6) 
^ '^^ * K2,i{x,y) K2,2{x,y) ' ^ ' 



with 



Ki,2{x,y) = i:"Ufi{x)M-J{eQ{y) 

K2,2ix,y) = -eix,y) + Y.lUiemx)M-l{ef)){y) 

provided that M is invertible, and {efi){x) = J-^ s{x,y)fi{y)diJ,{y). Note the 
order of indices in M~/ . 

4.2 Linear statistics 

Let a* and aj, j e be the creation and annihilation operator for the 
fermions and |0) be the state without fermions. The initial state is then 
given by 



i^^in)= n (4-8) 

j=-2N+l 

and the final state is 



l^fi-)= E n S*+2nJ0) (4.9) 
neCN j=-2N+l 

where Cn = {{no, ■ ■ ■ ,n-2N+i}\nj > rij-i^rij > 0}. Let us define the up- 
jump operator as 

ai = ^al^^ak, (4.10) 

which when applied on is actually a finite sum. Then the evolution from 
the initial state {t = 0) to the final one {t = Ty/2) is given by the transfer 
operator 

exp(Tai), f=^/2^T = 2T. (4.11) 
The linear statistics, i.e., for a bounded function g( : Z — > R, is 

E..r n (i-,(4^)))^i^MWi^^^^ (4.12) 
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where the x^"^, j G {— 2A^ + 1, . . . , 0} are the position of the fermions at time 
T^/2. Let us denote by p^"'\xi, . . . , x„) the n-point correlation function of 
Cjf"". Then 

\j=-2N+l ^ n>0 ' xi,...,x„e'Z j=l 

(4.13) 

For finite A^, p(") = forn > 2N. 

Proposition 4.1. Let us define the matrix $ with entries 

= T^—r:f^-'e{x - 2) , (4. 14) 

[x — ly. 

with <d the Heaviside function, the antisymmetric matrices S and A 

I + sgn{x - y){-lY I - sgn{x - y){-l)y 
Sx,y = ^ ^ sgn(x -y), (4.15) 

A, = J2 KAy%,j = [^'S<!>],,. (4.16) 

x,y€'E 

Then the n-point correlation function, for n G {0, . . . , 2N}, are given by 

p^^\x,, . . . ,x„) = Pf [K{x,,x,\^^,_^ (4.17) 

, . , • , , r^/ \ f Kii(x,y) Ki2(x,y) \ .,, 

where K ts a 2x 2 matrix kernel, K[x,y) = \ , i / ( , with 

Ki,2{^,y) = 'El=-2N+iK^Ai:}[^'s%,y = -K2M^) 

K2,i{x,y) = ~Y.l,=-2N+i\^'S%,^AT^^^^ 

K2,2{x,y) = Si^y-El=-2N+i[^'S\.A-j[<l>'S\y. 

When — > oo, ()4.13p becomes a Fredholm Pfaffian, Pf (J — Kg) = 

^^/Det{l — J^^Kg), where J = ^ ^ ^ , see Section 8 of |16^ . In this 

case, we consider bounded functions g with support bounded from below, so 
that the sum in ()4.13|) is well defined. From the point of view of operators, the 
determinant has to be though as defined through the modified determinant 
like in the case of the GOE case, see discussion at the end of Section 12.21 
Finally, note that A is invertible because Det(74) is the partition function of 
the line ensemble. 
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Proof. Since it is oi^en used, we denote the ordered set / = {— 2A^+ 1, . . . , 0}, 
and instead of writing a matrix M = [Mj ,,]jj=_2Ar+i,...^o we write M = 
[Mij]ij(zj. Let {x^'}), n G Cat as given in (jiT9|l . be the weight 

of fermions starting from positions {x™} = = i, and ending at 

{x^} = {x^^)j^i, = j + 2nj. The non-intersection constraint imphes [Ej 
that the weight can be expressed via determinants, 

HW''} - {4°}) = Det[^,,],,,, (4.19) 

with 

= (0|a,+2n,e^"^a*|0) = [e^"^]^.^^,^.,, = $,+2n,.. (4.20) 

Taking into account the even/odd initial position of the fermions, ()4.19|) can 
be rewritten as 



w{{x'^}-^{xl^}) = Det[Mxf)h^i n ^i4>i4U) (4-21) 

j=~N+l 

with 

e{x) = \ . o{x) = (4.22) 

Let us denote by p{x-2N+i, ■ ■ ■ ,xq) the probabihty that the set of end 
points {x^^.j = —2N + 1, . . . , 0} coincide with the set {x_2Ar+i, . . . , Xq}. We 
want to show that this probabihty can be written as a determinant times a 
Pfaffian. Since the x/s do not have to be ordered, let vr be the permutation 
of {—2N + 1, ... ,0} such that < X7r(i+i), that is, X7r(i) = xf'^, i & I. 
Moreover, define the matrix S = by setting Sjj = 5i,7r(j)- Then 

[M4nhei = [Mxj)h^i^- (4.23) 

Now let us show that 



n e(xf^)o(xg_i) = Pf (4.24) 

j=-N+l 

Since xf^ < xf^i, the components i,j {i < j) of the r.h.s. matrix are given 
by o(xf°)e(x5''). The Pfaffian of a matrix M = [Mij]ij(.i is 



Pf(M)= (-1)''^' n ^'^.^-i.-..' (4-25) 

o- i=-N+l 
0"2j-l<0"2i 

where the sum is on the permutations a of {—2N + 1, . . . , 0} with o"2i-i < 
(T2i. The identity permutation gives already l.h.s. of ()4.24|1 . Thus we have 
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to show that all other terms cancels pairwise. Take a permutation a such 
that a{2i — 1) < cr(2j — 1) < a{2i) < cr{2j) and define the permutation 
cr' by setting cr'{2j) = cr(2i), cr'(2i) = a'{2j), and <j'{k) = a{k) otherwise. 
The term of the Pfaffian coming from cr and a' are identical up to a minus 
sign because (—1)'°"' = —(—1)''^ '. Moreover, the only permutation for which 
a{2i - 1) < (t(2j — 1) < a(2i) < <j{2j) can not be satisfied for some i,j is 
the identity. Consequently ()4.24|) holds. 

Finally, define the matrix G = S* [S^.. Replacing the definition 

of S we obtain G = [S'*fl„ ^fin]i,jei- Then 

p(x_2jv+i, . . . , xo) = w{{x'^} ^ {xl^}) (4.26) 



= Det[$.(xf)],,,, n 

j=-N+l 

= Det[$,(a:,)],,,e/Det(S)Pf(S* [S^^^h^iE) 
= Det[$.(x,)],,,.7Pf[<,.J,,e/ 

where we used the property of Pfaffians Pf(E!*rS) = Pf(T)Det(S), see 
e.g. [22], and Det(S) = (-1)1^1. 

The probability ()4.26j) is of the form ()4.2|1 with 

e{x,y) = Sly, /.(a;)=<l>,,i (4.27) 

from which follows that 

M,,, = -A^„ {ef,){x) = -[S<!>l,i, (4.28) 
and the kernel is given by 

K\x,y)=( -l^^'f^y^ ^^'^^f'^W (4.29) 

But K and K' are two equivalent kernels (they give the same correlation 

functions) since K' = U^KU with U = i (^^^ ^ and Pf[U^KU] = 

Det[?7] Pf [i^]. We use K instead of K' uniquely because another derivation 
of the kernel gave K and we already carried out the analysis. □ 

5 Kernel for finite T 

In this section we compute the components of the kernel given in ()4.18|1 . At 
this stage we take the limit —>■ oo. The justification of this limit is in the 
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end of this section. The first step is to find the inverse of the matrix A. First 
we extend A to be defined for all j G Z by using ()4.16p to all Let us 
divide = where = {1, 2, . . .} and Z_ = {0, -1, . . .}. 

The inverse of A in ()4.18j) is the one in the subspace £^(Z_). Let us denote 
by P- the projector on Z_ and P+ the one on 'Z*^_. 

Lemma 5.1. The inverse of A in subspace P'i^J), which can be expressed 
as P^{P^AP^ + P+)"^P_, is given by 

[A-% = [a_ie-^"-P_e-^"^ - e-^"-P_e-^°^ai],,, (5.1) 
where = ^ij+i and «_i = a\. 

Proof. First we rewrite A as a sum of a Toeplitz matrix plus the remainder. 
Let be the matrix with [ae\i,j = ^ij^^^^'j^ ^ind = 1 — Og. Then 

S = Y, ^'«o - Yl "-'i^'"- (5-2) 

fc>0 fc>0 

It is then easy to see that, for Ve{x) an even polynomial of arbitrarily high 
order 

Ve{a±i)ae = aeVe{a±i), Ve{a±i)ao = aoK(a±i) (5.3) 
and for V^(x) an odd polynomial of arbitrarily high order 

Vo{a±i)ae = aoVo{a±i), Vo{a±i)ao = aeVo{a±i). (5.4) 

Hence A can be written as 

A = exp(Ta_i) ^(af - a^\^^ae){cosh{fai) + sinh(Tai)). (5.5) 

fc>0 

We pull the last factor in ()5.5p in front of the sum using the commutation 
relations ()5.3j) and ()5.4j) . and, after some algebraic manipulations, we obtain 

A = M + R (5.6) 

where M = i$*(g - g*)<l>, P = i(g + g*)(a„ - ae), with Q = Efc>o«f 
and $ = exp(Tai). 

Let B = [^~'^Y{a^iP^ - P_«i)$"^ We want to prove that it is the 
inverse of A in the subspace £^(Z_). First notice that Bij = if i > 1 or 
j > 1, which implies [A ■ B]ij = [P_AP_ ■ B]ij for i,j < 0. Therefore, for 
^,J<0, 

[A ■ P],„, = [(M + P) ■ [$-i]*f/o<l>-^].,, (5.7) 
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with 

Uq = q:_iP_ — P-tti, 
and, expanding M + R, we have 



-Tai 



-Tai 



+ 



-Tai 



(5.8) 



(5.9) 



where L^i = ^{Q — Q^)Uq and t/2 = \{.Q + Q*)(«o — «e)f^o- The components 
of these matrices are given by 



[U2] 



5„,ml[n<0] + ^5m,0 Sgn(n - 1)^^— ^-y^ 
1 l + (-l)" 



(5.10) 



2 ' 2 

and a simple algebraic computation leads then to [A ■ B]ij = 6ij for i,j < 0. 
Finally, since A and B are antisymmetric, [i? ■ A]ij = [A* ■ -B*] j,i = [A ■ B]j^i = 
6ij too. Therefore B is the inverse of A in the subspace £^(Z_). □ 

The second step is to find an explicit expression for the kernel's elements. 
Using the fact that of Lemma (5.11 is zero for i > 1 or j > 1, we can 

extend the sum over all i, j G Z and obtain 



x,y, 



-K2,iiy,x), 



Put ^ = e^'"ie-^°-^ We write S as in ()5.2|1 . use the commutation relations 



K2,iix,y) 

K2,2ix,y) = Sly-[S<^A-WS%,y. 

()5.3p and ()5.4|) . and after some straightforward algebra obtain 

K2,l = -^*(5f/o - f/i)^* - ^f/i^*, 
K2,2 = S' + SKl^, 

where Ui is given by ()5.1()|1 . and 

[^o]„,rrt ~ ((^n.m-l ~ '^m,n-l)l[n,m<0] 
1 1 + (-1)" 



(5.11) 



(5.12) 



[SUo - U^] 



m,0 • 



(5.13) 



2 2 

Using these relations we obtain the kernel elements, which are summed 
up in the following 
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Lemma 5.2 

with 



K{x,y) = Gix,y) + R{x,y), (5.14) 
Ri,iix,y) = 0, 

(-i)y 

Ri,2{x,y) = -^^J,+i(2T), 
R2,i{x,y) = ^^J,+i(2T), 
R2,2{x,y) = -S{x,y) + ^sgn{x - y) 

J2 Jy+2m{^f) + ^2n+x(2T), (5.15) 



2 ^ '^^^'"^ ' 2 

rra>l n>l 



Gi,i(x,i/) = -J2 Jx+n+i(2r)J,+„(2T) + ^ J,+„+i(2r)J,+„(2T), (5.16) 



n>l n>l 



Gl,2{x,y) = Jx+n{2f)Jy+n{2f) - J.+l(2T) ( J2 Jy+2m-l{2f) - M , 

n>l ^m>l ^ 

(5.17) 

G'2,i(x,i/) = - 5^ J,+„(2T) J,+„(2r) + J,+i(2r) ( Jx+2m-i{2f) - ^ j , 

n>l ^m>l ^ 

(5.18) 

</x'+2m(2T) Jj^+2n+l(2T) — J^+2m+l (2T) Jjy+2n (2T) 

m>l n>m ?i>l m>n 

Yl J-+^rn{2f) + V2n(2T) " \ sgn(x - y), (5.19) 



m>l 7i>l 



where Jm{t) denotes the rrfi^ order Bessel function. 

Remark: this result could also be deduced starting from Section 5 of [12] . 
Now we justify the N oo limit. Let us first explain the idea. Denote 
the sets / = {-2N + 1, . . . , 0} and L = {-A^ + 1, . . . , 0}. We consider 
the kernel's elements for x,y > 0. For {i,j) G \ L^, the inverse of A 
for finite N differs from the inverse for N = oo only by 0{e~^^) with fi = 
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fi{T) > 0. On the other hand, the contribution to K^ {x,y) coming from 
€ are exponentially small in A^. Therefore, replacing the inverse 

of A for finite N with the inverse obtained in Lemma f5. II we introduce only 
an error exponentially small in N. The dependence of the kernel's elements 
on N is only via the extension of the sums in ()4.18j) . which limit the one we 
derived in Lemma f5. 21 

In what follows we denote hj Aj^ the 2N x 2N matrix ()4.16|) and by A 
the N = oo one. 

Lemma 5.3. // we replace [Aj/']ij by A~j in the kernel's elements d^. j<§| ), 
then for N large enough, the error made is 0{e^^^) for some constant fi = 
fi{T) > 0. The error is uniform for x,y > 0. 

Proof. Here we use some results of Appendix lA.il First, we define the matrix 
B by setting, Bij = Arj for (i, j) e I x L, and Bij = _2N+i-j for 

e / X (/ \ L). Since [A^jij = -[A^] by jQ) follows 

that 

AnB = 1-C (5.20) 

for some matrix C with ||C|| = maxjj \Ci,j\ < 0{e~'^^'^). Therefore, for A^ 
large enough, 

A],^ = B{1 + D), D = ^C^ (5.21) 

fc>i 

with \\D\\ < (9(e~^2^) too. Thus, replacing A^^ with B we introduce an 
error in the kernel's elements of 0{N'^e~^'^^). 

If we replace Bi^j with A^^ also in (i, j) E L x (I \ L) we introduce an 
error of (9(A^^e~^^^), with /i3 = min{/ii, /i2/2}. This is achieved using ()A.6|1 
for i < j + N/2, and ()A.4jl otherwise. 

The final step is to show, using only the antisymmetry of Ajj^ that the 
contribution of K coming from G (/\ L)^ are also exponentially small 
in A^. For (i, j) G (/ \ L)^, it is easy to see that, uniformly in x,y > 0, 

= 0(e-^^^) (5.22) 

(5$)^,i = 0(e-^^^), for odd X, 

{S(^)x,i = sinh(T) + C(e^'^^^), for even x and even i, 

{S<^)x,i = cosh(r) + O(e-^i^), for even X and odd 2. 

Therefore, the contributions for Ki i, K12, and K2^i are 0{N'^e~^^'^) be- 
cause they contain at least a factor (9(e~'^^^) coming from $^ j or For 

i^2,2 there are terms without (9(e~'^^^), and containing only sinh(T) and/or 
cosh(T). These terms cancel exactly because Aj^^ is antisymmetric. Conse- 
quently, we can simply replace Bij with A~j also in (i, j) G (/ \ L)^ up to 
an error 0{N'^e-''^^). ' □ 
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6 Edge scaling and asymptotics of the kernel 



In this section we define the edge scahng of the kernel, provide some bounds 
on them which will be used in the proofs of Section [3 and compute their 
T ^ oo limit. 

The edge scaling of the kernel is defined by 

^$'(6,6) = f'/'G,{[2f + ^,f'/%[2f + ^,f fc=(l,2),(2,l) 
GT?2(ei,e2) = G2,2i[2f + ^^f'/%[2f + ^2f'/']), (6.1) 

and similarly for R'^^{^i,^2)- 

Next we compute some bounds on the kernel's elements such that, when 
possible, they are rapidly decreasing for i^i, 1^2 ^ 1- 

Lemma 6.1. Write 

X <0 Q ( \ _ i + kl, a; < . . 
^^o[x) - I exp(-x/2), a; > ' '^^^ ~ \ exp(-x/2), x > ' 

^^(-) = l ^'V^'^o. "^n • (6-3) 
' \ exp(-x/2), X > ^ ' 

Then there is a positive constant C such that for large T 

i<f:2(ei,6)i < ci]o(6), 

i<i:i(ei,6)i < cno{i2). (6.4) 

|i?T;i:2(ei,6)l < C(l]i(6)+f^l(6)), 

and 

|Gg:,(ei,6)i < ^7^2(6)^^2(6), 

|Ggyei,6)i < cfii(6)(i + f^2(6)), (6.5) 

|Gg:,(ei,e2)l < c^]l(6)(l + ^^2(6)), 

|G^T?2(6,6)l < C(i + fii(6) + ^^i(6) + fii(6)fii(6))- 

Proof. We use Lemma fA. II and Lemma fA. 21 to obtain the above estimate. 
1) The bounds on |_R^f^2(65 6)1 and |i?Ij?|'^]^(6, 6)1 are implied by 
Lemma lA.ll 
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2) Bound on li^^^f^^i, ^s)!- 

I^T?2(ei,6)| <\ + \ J2 l^[2f+(2M+6)TV3](2T)| + ^ 6) (6.6) 



4 2 



and 



(2T)|. (6.7) 

MeN/Ti/3 MeN/fi/3 

For ^2 < 0, 

(EH) < ^ l'^[2f+Mri/3](2T')| + ^ I J[2j.+j,,jfi/3](2T)|. 

Me6+IN/Ti/3uK2,0] MeN/fi/3 

(6.8) 

By ()A.8j) the first term is bounded by a constant times (1 + |^2|) and by ()A.9|) 
the second term by a constant. For ^2 > 0, 

m< E l'^[2T+MTV3](2f)| (6.9) 

Mea+N/Ti/3u[^2,oo) 

which, by ()A.9j) . is bounded by a constant times exp(— ^2/2). Therefore 

E I V+(Af+6)TV3](2T)| < Cn^i^2). (6.10) 

MGlN/f 1/3 

for a constant C, from which follows the desired bound. 

3) Bound on |G^f^^(^i, ^2)!- Let us define J„(t) = J„+i(t) — J„(t). Then 

^T;l,l(^l' '^2) = ^"^^^ E '^[2T+(a+M)Ti/3](2T) Jj2'f+(^2+M)Ti/3](2T) 

Af6N/fi/3 

- (6-- 6)- (6.11) 

For large T, the sums are very close integrals and this time we use both 
Lemma fA. II and Lemma fA.2t obtaining 

POO 

|GT'fi(ei,6)| < C dMno{M + ^^)n,{M + ^2) 

Jo 

POO 

< C dMni{M + ^i)ni{M + ^2) (6.12) 



for a constant C > 0. It is then easy to see that r.h.s. of ()6.12|1 is bounded 
as follows: for 6 < 6 < by C(l + \^^\)\ for ^1 < < ^ by C(l + 



25 



|ei|)'exp(-6/2), and for < < 6 by Cexp(-ei/2) exp(-G/2), for some 
other constant C > 0. Therefore IGr^il^i^ 6)1 < ^^^^2 (6)^2 (6)- 

4) Boundon|Ggyei,e2)|. 

^tIi^^'^) = '^[2T+(6+M)fV3](2T)(r^/^Jj2f+(^2+M)Ti/3](2T)) 
AfeN/Ti/3 

•^[2f+{Ci+M)f 1/3+1] (2T) ^ ^ '^[2T+(6+2M)fi/3_i](2T) — - 

AfeN/fi/3 

In the first sum, the term with ^2 is bounded by a constant and remaining sum 
was already estimated in ()6.10|1 . The second term is bounded by a constant 
times f2o(^i)f^2(6)- Using fio(6) < ^1(6) "we conclude that \G'^^\{^i,^2)\ < 

cni(ei)(i + fi2(6)). 

5) Bound on ^2)1- The bound is the same as for \G'j..l\{^i,^2)\- 

6) Bound on \G'^.2\i^i^^2)\- The terms with the double sums are esti- 
mated applying twice ()6.1Up and are then bounded by ^1(^1)^1 (^2)- The 
two terms with only one sum are bounded by fii(G) and fii(.C2) respec- 
tively, and the signum function by 1/4. Therefore, for some constant C > 0, 

IG^'iy^' 6)1 < c(i + ^]l(el) + ^^i(6) + (6)^^1(6))- □ 

Finally we compute the pointwise limits of the G"s since they remains in 
the weak convergence. 

Lemma 6.2. For any fixed ^1, ^2; 

lim G?;r(ei,6) = ^^"'^(ei,^), (6.13) 

where the 's are the ones in 112. 

Proof. Let us consider ^1,^2 fixed. In the proof of Lemmas 16.11 we have 
already obtained uniform bounds in T for G'^|'^(.^i, ^2)5 so that dominated 
convergence applies. To obtain the limits we use ()A.33|) . i.e., 

^lim Ti/V[2T+5TV3](2T) = Ai(0, (6.14) 

and 



^hm T2/3(Jp^+^^,/3+i](2T) - J[2T+€Ti/3](2T)) = Ai'(e). (6.15) 

The limit of G^'|?'i(6,6) follows from (IFTTTD . 
The limit of ^^^^2(6? 6) leads to 

/•OO 1 / /"OO \ 

dAAi(ei + A)Ai(6 + A)--Ai(ei)( y dAAi(e2 + A)-lj (6.16) 
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which equals Gf^^^ since J^^ dX Ai{^2 + X) - I = - dX Ai(6 - A). 
The hmit of G'^.2\{^i,C,2) is obtained identically. 
Finally, the limit of G^l^jl^b^a) is given by 

-| roc roo 

- dX d/i Ai(ei + A) Ai(6 + /x) - (6 ^ 6) (6.17) 
dAAi(ei + A) + - / d/iAi(e2 + /i)-7Sgn(ei-6), 



which can be written in a more compact form. Since J^dA Ai(A) = 1, 

/ dAAi(ei + A)=/ dA / d^Ai(ei + A)Ai(e2 + /i), (6.18) 

Jo Jo J~oo 

and the signum can be expressed as an integral of Ai(^i + A) Ai(^2 + f^) 

-sgn(ei-e2)= / dA / d/iAi(ei + A)Ai(e2 + /i)sgn(A-/x). (6.19) 
Jr Jn 

In fact 

r.h.s. of (j6T9|l = / dA / d/iAi(A)Ai(/i)sgn(A -/i + C) = KC) (6.20) 
Jn Ju 

with C = ^2 ~ ^1- For C = it is zero by symmetry. Then consider C > 0, the 
case C < follows by symmetry. By completeness of the Airy functions, 



[ d/iAi(/i)Ai(/i-C)=5(C). (6.21) 
Jn 



dC 

Then using ()6.18|) and ()6.19p we obtain the result. □ 



Remark that the GOE kernel in differs slightly from the one written 
here, but they are equivalent in the sense that they give the same correlation 
functions. 

For the residual terms the limit does not exist, but exists in the even/odd 
positions. In particular 

1 /"°° 

E '^[2T+«tvh2h(2^) = 2 / '^^^'^^ + (^-22) 
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7 Proof of Theorem 12.11 

In this section we first prove the weak convergence of the edge rescaled point 
process of r^^™ to 7^*^oe ^^^q y ^ qq hmit. Secondly, using the equivalence 
of the point process and Ct^^^ '^6 prove Theorem 12.11 

Theorem 7.1. Let us define the rescaled point process 

vTU) = E /((^ - 2T)/f V3)c^y-(a;) (7.1) 

with T = ■\/2gT = 2T and f a smooth test function of compact support. In 
the limit T ^ oo it converges weakly to the GOE point process, i.e., for all 
m G N, and fi, ■ ■ ■ , fm smooth test functions of compact support. 



(7.2) 



where the GOE kernel is given in 



Proof. Let fi, . . . , fm be smooth test functions of compact support and 
/i(x) = /.((x-2r)/Ti/3)^then 

Erf ]^?7T™(/fc) ) = fl{xi) . . . fm{Xm)Pi[K{Xi,Xj)]ij=i^,„^rn 

^ k=l ^ xi,...,Xme'Z 

— T Y fi{xi) . . . fmiXrn)Pf[L{Xi,Xj)]ij 

=l,...,m 

(7.3) 



Xi,...,XmeIj 



f T^/^ \ 

where X = \ ^ ^1 and L{x,y) = {XKX^){x,y), i.e., Li^i{x,y) = 

f2/3Ki,i(x,y), Lkix,y) = fi/3irfc(x,y), for k = (1,2), (2,1), and L2,2ix,y) = 
K2^2{x,y)- Moreover, we define the edge scaling for the kernel elements as 

^Tlr(6,e2) = L,{[2f + ^,f'/% [2T + 6T^/3]). (7.4) 

In what follows we denote by = (xj — 2T)/T^^^. To simplify the notations 
we consider T e N, but the same proof can be carried out without this 
condition, replacing for example Z/T^/^ (^ _ 2T)/T^/3 [y^ (j7.5j) . Then 



(7.5) 
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Let us denote = [^iT^/^]/T^^^ the "integer" discretization of ^j. Then 

Et ( n "(A)) = /^^ d6 ■ ■ ■ dai(eO ■ ■ ■ uo ^iiL'^^'^iU ej)].,=i,....™. 

(7.6) 

Using the definition in (jfj.lj) we have 

L^r (6, 6) = G?/ (ei, 6) + i?Tlr (6, 6), (7.7) 



therefore ()7.(i|l consists in one term with only G^f^ plus other terms which 



contain at least one R'^r^!^'^ 



First consider the contribution where only Gf^.f^ occur. Let Mj > be 
the smallest number such that fj{x) = if > Mf, for all j = 1, . . . ,m. 
We bound the product of the /j's by 

m 

i/i(e()---/™(e^)i <nii/^-ii-i[-M/>^/i(^^-) (^-8) 

and, in the same way as in Lemma 17^ but with K^f replaced by G^^, we 
conclude that this is uniformly integrable in T. We then apply dominated 
convergence and take the limit inside the integral obtaining 

= [ d^i ■ ■ ■ dUH^l) ■ ■ ■ fm{U) Pf [G^''''(ei, Q]i,3=l,...m. (7.9) 

Next we have to show that whenever some R'^r^^ are present their contri- 
bution vanish in the limit T oo. In ()7.(i|) we have to compute the Pfaffian 
of Et defined by 



^ L^^f'i((n + l)/2,(/ + l)/2), n odd, /odd, 
^T^yi^ + l)/2, V2), n odd, / even 



L^^f,(n/2, (/ + l)/2), n even, / odd, ^ '''^^^ 

L^|'^2(^/2, //2), n even, / even. 



for 1 < n < / < 2m, with L^^f (a,6) = Lf^^{Uib)- The Pfaffian of Et is 
given by 

Pf(ET)= (-l)'"'^T(al,^T2)■■■ET(^X2„^-l,a2„,). (7.11) 
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Now we have to check that the product of residual terms does not contain 
twice the term (—1)^ for the same x. This is imphed by Lemma (7. 2[ 

Let us decompose the sum in ()7.5|) into 2™ sums, depending on whether 
is even or odd. Denote ^| = [^if^/^ /2]2/f^/^ and = {[^if^^^ /2]2 + 
1)/T^/^ the "even" and "odd" discretizations of ^j. Then 

m = 7^ J2 I d6-- -dewier) ■■■^(e)pf[^T^^(c^e;oL=^^^^^^^^ 

Si={o,e], 
i=l,...,m 

(7.12) 

With this subdivision, each term in the Pfaffian converges pointwise to a well 
defined limit. Moreover all the 2*" integrals, including G^f'^'s and/or i^^'^l'^'s, 
are uniformly bounded in T. By dominated convergence we can take the 
limit inside the integrals. 

Each time that there is a R'^.l\{^i, , or R^^\{^j,^i), the integral with 
Si = o and the one with Si = e only differs by sign, therefore they cancel each 
other. Each time that appears RT^^\{C,i,C,j), the part including coming from 
the (—1)^* and the one with (— 1)^^ simplifies in the same way. Finally we 
consider the second part, the one including the S and signum function. The 
sum of Si = {o,e} and Sj = {o,e} of the terms with —S{^iT~^^^,^jT~^^^) 
equals minus the ones with isgn((^j — ^j)T^^/^). Consequently all the terms 
including at least one time -R^f' have a contribution which vanishes in the 
T — i> oo limit. □ 

Lemma 7.2. The following products do not appear in \7.11\ ): 

Proof. We prove it by reduction ab absurdum. We assume that the product 
appear and we obtain a contradiction, (a) appears if there exist some a < b 
and c < d with a, b, d even and c odd, all different, such that i = a/2, j = 6/2, 
k = (c+ l)/2, i = d/2. But this is not possible since d ^ a. (b) appears if 
there exist some a < b and c < d with a, b, d even and c odd, all different, 
such that i = a/2, j = 6/2, k = {c + l)/2, j = d/2. But this is not possible 
since d ^ b. (c) appears if there exist some a < b and c < d with a, b, c even 
and d odd, all different, such that i = a/2, j = 6/2, i = c/2, k = {d + l)/2. 
But this is not possible since c ^ a. (d) appears if there exist some a < b and 
c < d with a,b,c even and d odd, all different, such that i = a/2, j = 6/2, 
j = c/2, k = {d + l)/2. But this is not possible since c ^ b. (e) appears if 
there exist some a < b and c < d with 6, c even and a, d odd, all different. 
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such that i = {a+ l)/2, j = b/2, j = c/2, k = {d + l)/2. But this is not 
possible since c 7^ 6. □ 

Lemma 7.3. There exists a constant C > such that 

lET{\vTiM-M,oo))r) < C^e^'^l\mTl'' (7.14) 

uniformly in T. 

Proof. The m-point correlation function p^™^(.Ci, . . . , ^m) is a sum of product 
of K^l'^^s which contains twice every ^j's, i = 1, . . . , m, and only in K!^f^ 
the two argument can be the same. From Lemma f6. If for any ^1,^2 ^ IR, 

mti(^^,^2)\ < Cexp(-6/2)exp(-6/2) (7.15) 
|<f:2(ei,6)l < C^exp(-6/2) 
\K^r%(^„^,)\ < Cexp(-6/2) 

I^Ste,^)! < c. 

For negative ^ we could replace exp(— ^i/2) by (1 + |^i|)^ where appears, but 
for our purpose this is not needed. 

All the products in p*-'"^(^i, . . . ,^m) contain at least one exp(— ^j/2) for 
each i. In fact, this holds if: K!^2\{^1j^2) is not multiplied by i^^^'^f 2(^35 ^2), 
K^%i^,,^,), <i:i(ei,e3), <f:i(6,e3), and if <f:2(6,6) is not multi- 
plied by -ft'^tl'il^a, ^3)- This is already proven in Lemma 17^ 

Consequently, 

ET(ir/r(i[-M,oo))r) = / dei...de™p?")(6,...,u 

^ ^ J[-J\/,oo)'" 

< {2mr/^( [ Cexp(-^/2)d^)"' = 2™C""e*^'"/2(2m)"^/2(7.16) 

^ J[-M,oo) ' 

uniformly in T. The term (2m)™'/^ comes from the fact that the absolute 
value of a determinant of a n x n matrix with entries of absolute value not 
exceeding 1 is bounded by n"/^ (Hadamard bound). Finally resetting the 
constant as C2\/2 the lemma is proved. □ 

To prove Theorem 12.11 we use Theorem 17. H Proposition 13. 2[ and 
Lemma 17.31 

Proof of Theorem \2.1\ Let us denote by Xj, j < 0, the position of the j**^ 
element of C,^^^ and a;^'*'™, j < 0, the position of the j'*^ element of Ct^- Then 
define ^j^T and by 

Xj = 2T + ij^TT^'^2-^'\ xf"" = 4T + (2T)i/l (7.17) 

31 



By Proposition 13.21 Xj — j = ^ j); which imphes 



J 



Let fi,...,fm be test functions of compact support and denote by Mf > 
the minimal value such that fj{x) = if |x| > Mf, j = 1, . . . ,m. Then 

m m 

ET(n^T"(A)) = ET(nE^(C" + V(2T)^/^)) (7.19) 

fe=l k=l i<0 

m 

= E n/^(c,T+w(2T)^/=^)). 

ii,...,im<0 fc=l 

We bound the f^s by their supremum times 1[-M/,M/] as in ()7.8|1 . then 



Ir.h.s. of (mnH < E^(n E lM//.M,](e-r + V(2T)^/^)) n ll/.lloo, 

(7.20) 

and, since IL[_m„m,](CT + ^(2^)1/^) < l[_M„oo)(err + V(2T)V3) < 
l[-A'//,oo)(C\'r')5 it follows that 

m m 

Ir.h.s. of (Eisiii < ET(n^T"(i[-A/„oo))) n 11/^11- c^-^i) 

i=i i=i 

which is uniformly bounded in T from Lemma 17.31 Therefore by Fubini's 
theorem, 

m m 

ET(n^T"(/.))= E E^(n/^(CT + W(2T)i/^)). (7.22) 

fc = l ii,...,jm<0 fc = l 

Moreover, UifJ^ + 2./(2T)V3) = + /^(e,„^)^,/(2r) V3 for some 

4,T e K^,T + ^/^/(2T)^/^e:^,T]• Therefore dm equals 

m m 

E ET(n^(c,?))=ET(n^r(A)) (7-23) 

il,...,im<0 fc=l fe=l 

plus 2™ — 1 terms which contains some /^(^jj.^T)^fc/(2T)^/3. Finally we have 
to show that these terms vanish as T — > oo. First we bound the /^'s and 
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the /^'s by ||/fc||oo and ||/^||oo times 1[_m^,M/]- Therefore each of the 2™ — 1 
terms is bounded by a 

^ mm 

^nii/^iunii/^ii- E ni^^iE^(niw-)(c^,T)) (7.24) 

kel k€J ii,...,im<0 keJ k=l 

where I and J are subset of {1, . . . , m} with / U J = {1, . . . , m} and J is 
non-empty. Let jo = min{zi, . . . , im}, then 



ET(ni[-M/.oo)(C,T)) = Er(lhM„oo)(C',T); 
P^(C't > -Mf) < PT{vriM-M„oo)) > Jo) (7.25) 



since |jo| > for all + 1, . . . ,m. From ()7.25j) it follows that ()7.24j) is 
uniformly bounded in T and vanishes as T — oo. We have then proved that, 
for all fi, . . . , fm smooth test functions of compact support, 

m mm 

^hm E^(n4^*(/.)) =^li^E^(nr/r(A)) =E(n^^°V;.)), 

k=l k=l k=l 

(7.26) 

the last equality being Theorem 17.11 
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A Appendices 

A.l Bounds on the inverse of A 

Let us denote the finite matrix Ahj and its inverse by Aj/^. For the N = 
oo case we use the notations A and A~'^. Let us denote / = {— 2A^ + 1, . . . , 0} 
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and L = {—N + 1, . . . , 0}. Using ()5.6|) we have 

r-pk—i ^l—j ]^ 



To obtain some properties of A ^ , we first estimate [e ^°'-^P_e '^"^jij 

{k-t)\ {k-jy. 



max{i j'}<A;<0 



(-l)l*-^l/|,_,l(2T)- J] ^ ^ 



where Jfc is the modified Bessel function / of order k. From ()A.2jl and (Z + 

\i-j\V- > - j|! follows 

|[e-^"-P_e-^-],,| < /o(2T)^- < ^-e^^ (A.3) 
which implies 

\Arj\ < 2^^e'^ < ci(T)e-'^^(^)l-^-| (A.4) 

for some constants Ci,/i2 > 0. 

The remainder sum in ()A.2j) is exponentially small in — max{i,j}. In 
fact, for n = — max{z,j}, 

|[e-^"-P_e-^-].,-(-l)l^-^-|j|._,l(2T)| 

for some constant fii > 0. Thus, for all (i, j) such that max{i, j} < —N/2, 
\Arj - hm < C2(T)e-'^^^/2 (A.6) 

for some constant C2 > 0, that is, in this region A~^ is exponentially close to 
a Toepliz matrix. 

For j G L, using (jA.lj) and ()A.4|1 . we obtain 

\{A^A~^ - = A^iAll < C3{f)e-^^^ (A.7) 

1<-2N 

with C3 > a constant. 
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A. 2 Some bounds 
Lemma A.l. For N >0, 

|Ti/V[2T+jvri/3](2T)| < exp(-iV/2)0(l) (A.8) 
uniformly in T > Tq for some constant Tq. 

For N < it follows from a result of Landau Ol) /. see \A.34 ), that 

\T^'^J[2T+ntH'^T)\<C (A.9) 
uniformly in T for a constant C > 0. 

Proof. To obtain the bound we use 9.3.35 of P, i.e., for z G [0, 1], 



Jn{nz) 
where 



1-z 



^1/3 ^5/3 



(AlO) 



■I 2/3 



C(;z) = (3/2)2/Mn(l + v^r^^) - ln(2) - v^r^^ . (A.ll) 

In our case, n = 2T + NT^/^ and z = (1 + e)"^ with e = ^AT-^/^ > 0. 
This imphes that z G [0, 1]. In this interval the function (^{z) is positive and 
decreasing. The prefactor is estimated using A({z{e)){l — z{e)'^)~^2~^/^ < 
1 + |e for all e > 0. Moreover, for a; > 0, Ai(x) < Ai(x/2) and | Ai'(x)| < 
Ai(x/2). Therefore 

|Ti/VpT+^^V3](2T)| < (1 + Ai(n2/3C/2)(1 + OiT''/')) (A.12) 

where we also used {2TY^^ < n^^^. Next we bound ()A.12jl separately for 
N < iT2/3 and N > ^T^/\ 

Case 1) < A < ^T^/\ In this case e <l and, for e G [0, 1/3], C{z{£)) > ^ 
holds. Replacing n by 2T in the Airy function we have an upper bound since 
it is a decreasing function, consequently 

V+^Ti/3](2T)| < 2 Ai(A2-^/3)(^ ^ 0(T-^/3)). (A.13) 
Finally it is easy to verify that 2Ai(A2-^/3) 

< exp(— A/2), and obtain the 

bound of the lemma. 

Case 2) A > ^T^/^. In this case £ > ^ and z{e) < |. In this interval 
C{z) > \{\n{8e)Y^^ from which follows 

l^'^'V+^TV3](2T)| < (AT-=^/3)^/^Ai (|(nln(4AT-2/3))2/3) Q^iy (a.14) 

35 



For a; > 0, Ai{x) < exp{-lx^/^), and N > \T'^I^ implies N = ANT'^/^ > 2. 
Consequently, 

|r''"V+iVTV3](2T)| < 7Vi/^exp(-ciT(l + iV/8))0(l) 

< exp(-cir) exp(-2c2T7V)7V^/^0(l) (A.15) 

with ci = ln(2)/3, C2 = Ci/16. For T > 10 and TV > 2, TVV^ exp(-C2T7V) < 
1, and exp(— C2T'7V) < exp{—N/2) for T large enough. These two last in- 
equalities imply 

|7^'/'V+iVTV3](2T)| <exp(-ciT)exp(-iV/2)0(l) (A.16) 

for T large enough, and the lemma is proved. □ 

Lemma A.2. For all N >0, 

Dt,n = \T^'\J[2T+NTV^+ipT) - J[2T+iVTi/3](2T))| < exp(-Ar/2)0(l) 

(A.17) 

uniformly inT >Tq for some constant Tq. 
For N < 0, there is a constant C > such that 

Dt,n<C{1 + \N\) (A.18) 

uniformly inT > 1. 

Proof First we consider > 0. Let N' ^ N + T''^/^, then we have to 
subtract ^[2r+Arri/3](2T') to J[2T+Ar'ri/3](2T). In term of £ = \NT~'^/^ the 
difference is 1/ {2T) . Let us define 

and 

/(^)=(^Ai[(2T)2/3p(£)]. (A.20) 

With these notations, 

Ji2T^NTH2T) = f{e) + j^,M[{2Tn{e)]0{T-') 

+ Ai'[(2r)2/3p(£)]0(r-^/3). (A.21) 
Now we bound Dt n as follows. 
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Case 1) Let us consider N G [0, |T^/^]. The second and the third terms are 
simply bounded by their absolute value. Then 



Dtn\ < 



fie+^)-f{e] 



The first term is bounded by 



2^2/3 

where 



1 



2T' 



1 



f{e + —)-f{e) <T2/3 sup , (A.23) 



xG[e,e+l/2T] 



2T' 



(A.24) 

We are considering the case of G [0, |T^/3], which corresponds to £ G 
[0, 1/4]. The functions g, g', and q-p' behave modestly in this interval. They 
satisfy 

G [1.22,1.26], |g'(x)| G [0.14, 0.17], (x)| G [1.3, 1.6] (A.25) 

for X G [0,1/4]. The Airy function and its derivative are bounded as in 
Lemma fA. II Therefore 

\DtM < exp(-Ar/2)(l + 0(r-2/3)). (A.26) 

Case 2) Let us consider > \T^/^. This case is simpler. We apply (IA.16|1 
and obtain the bound 

\Dt,n\ < T'^'^ exp(-cir) exp(-Ar/2)C»(l) < exp(-A^/2)0(l) (A.27) 

for T large enough. 

Secondly we consider A^ < 0. For |A^| > T^/^, using ()A.34jl we obtain 

\DtA < csT^/^ < CslA^I (A.28) 

for some constant C3 > 0. Next we consider |A^| < T^^^. Since A^ is negative, 
z > 1 and (jA.lOj) holds with ({z) given by 

r 1 2/3 

((2) = -(3/2)2/3 - 1 - arccos(l/2) . (A.29) 
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Recall that z = {1 + e)'^ and e = iA^T'^/s. |e| < iT"^/^ is very close to 
zero. The estimate follows the same outline as for the case 1) for positive N. 
Take T > 1, then e G [-|, 0] and 

G [1.25, 1.37], G [0.16, 0.25], |g(£y(e)| G [1.5, 3.1]. (A.30) 

The difference is that now the Airy function in not rapidly decreasing since 
pi^) < and its derivative is even increasing. We use some simple bounds: 
I Ai(x)| < 1 and | Ai'(x)| < 1 + |x| for all x, with the result 

\Dt,n\ < C4(l + |iV|)(l + 0(r-2/3)) (A.31) 

for a constant C4 > 0. □ 

A. 3 Some relations involving Bessel functions 

Here we give some relation on Bessel function P which are used in the work. 
Bessel functions J„ are defined via the generating function by 

exp {lz{t - 1/t)) = J^t'^Jkiz), {t ^ 0). (A.32) 

Then 

1. for n G N, J_„(z) = (-1)"J„(2), 

2. Mz) + 2Zk>iJ2k{z) = l, 

3. J^{z)+2Z,>^Ji{z) = h 

4. for n > 1, ET=oi-'^)'M^)J2n~k{z) + 2 ^^=1 Jk{z)J2n+k{z) = 0. 
Moreover the limit 

^hm ri/V[2T+«^i/3](2r) = Ai(0 (A.33) 

holds. A useful result of Landau [TT] is the following: 

I J„(x)| < c|x|-^/^ c = 0.785... for all n e Z. (A.34) 
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